Tracing tumor heterogeneity of pleomorphic carcinoma of the lung.

BACKGROUND: Pulmonary pleomorphic carcinoma (PPC) is an aggressive and highly heterogeneous non-small-cell lung carcinoma whose underlying biology is still poorly understood. METHODS: Forty-two tumor areas from 20 PPC patients were microdissected, including 39 primary tumors and 3 metastases and the histologically distinct components were subjected to whole exome sequencing (WES) separately. We further performed in silico analysis of microdissected bulk RNAseq and methylation data of 28 samples from 14 PPC patients. We validated our findings using immunohistochemistry. RESULTS: The epithelial and the sarcomatoid components of PPCs shared a large number of genomic alterations. Most mutations in cancer driver genes were clonal and truncal between the two components of PPCs suggesting a common ancestor. The high number of alterations in the RTK-RAS pathway suggests that it plays an important role in the evolution of PPC. The metastases morphologically and genetically resembled the epithelial or the sarcomatoid components of the tumor. The transcriptomic and epigenetic profiles of the sarcomatoid components of PPCs with matched squamous-like or adenocarcinoma-like components differed from each other and they shared more similarities to their matched epithelial components. NCAM1/CD56 was preferentially expressed in the sarcomatoid component of squamous-like PPCs, whereas CDH1/E-Cadherin expression was downregulated in the sarcomatoid component of most PPCs. CONCLUSION: LUAD-like PPCs are mainly driven by RTK-RAS signaling, whereas epithelial-mesenchymal transition programs as highlighted by increased NCAM1 and decreased CDH1 expression govern the epithelial-sarcomatoid transition between the clonally related tumor components. Several alterations in PPCs pinpoint therapeutic opportunities.

High expression of HOXA13 correlates with poorly differentiated hepatocellular carcinomas and modulates sorafenib response in in vitro models.

Hepatocellular carcinoma (HCC) represents the fifth and ninth cause of mortality among male and female cancer patients, respectively and typically arises on a background of a cirrhotic liver. HCC develops in a multi-step process, often encompassing chronic liver injury, steatosis and cirrhosis eventually leading to the malignant transformation of hepatocytes. Aberrant expression of the class I homeobox gene family (HOX), a group of genes crucial in embryogenesis, has been reported in a variety of malignancies including solid tumors. Among HOX genes, HOXA13 is most overexpressed in HCC and is known to be directly regulated by the long non-coding RNA HOTTIP. In this study, taking advantage of a tissue microarray containing 305 tissue specimens, we found that HOXA13 protein expression increased monotonically from normal liver to cirrhotic liver to HCC and that HOXA13-positive HCCs were preferentially poorly differentiated and had fewer E-cadherin-positive cells. In two independent cohorts, patients with HOXA13-positive HCC had worse overall survival than those with HOXA13-negative HCC. Using HOXA13 immunohistochemistry and HOTTIP RNA in situ hybridization on consecutive sections of 16 resected HCCs, we demonstrated that HOXA13 and HOTTIP were expressed in the same neoplastic hepatocyte populations. Stable overexpression of HOXA13 in liver cancer cell lines resulted in increased colony formation on soft agar and migration potential as well as reduced sensitivity to sorafenib in vitro. Our results provide compelling evidence of a role for HOXA13 in HCC development and highlight for the first time its ability to modulate response to sorafenib.

MYBL1 rearrangements and MYB amplification in breast adenoid cystic carcinomas lacking the MYB-NFIB fusion gene.

Breast adenoid cystic carcinoma (AdCC), a rare type of triple-negative breast cancer, has been shown to be driven by MYB pathway activation, most often underpinned by the MYB-NFIB fusion gene. Alternative genetic mechanisms, such as MYBL1 rearrangements, have been reported in MYB-NFIB-negative salivary gland AdCCs. Here we report on the molecular characterization by massively parallel sequencing of four breast AdCCs lacking the MYB-NFIB fusion gene. In two cases, we identified MYBL1 rearrangements (MYBL1-ACTN1 and MYBL1-NFIB), which were associated with MYBL1 overexpression. A third AdCC harboured a high-level MYB amplification, which resulted in MYB overexpression at the mRNA and protein levels. RNA-sequencing and whole-genome sequencing revealed no definite alternative driver in the fourth AdCC studied, despite high levels of MYB expression and the activation of pathways similar to those activated in MYB-NFIB-positive AdCCs. In this case, a deletion encompassing the last intron and part of exon 15 of MYB, including the binding site of ERG-1, a transcription factor that may downregulate MYB, and the exon 15 splice site, was detected. In conclusion, we demonstrate that MYBL1 rearrangements and MYB amplification probably constitute alternative genetic drivers of breast AdCCs, functioning through MYBL1 or MYB overexpression. These observations emphasize that breast AdCCs probably constitute a convergent phenotype, whereby activation of MYB and MYBL1 and their downstream targets can be driven by the MYB-NFIB fusion gene, MYBL1 rearrangements, MYB amplification, or other yet to be identified mechanisms. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Bi-allelic alterations in DNA repair genes underpin homologous recombination DNA repair defects in breast cancer.

Homologous recombination (HR) DNA repair-deficient (HRD) breast cancers have been shown to be sensitive to DNA repair targeted therapies. Burgeoning evidence suggests that sporadic breast cancers, lacking germline BRCA1/BRCA2 mutations, may also be HRD. We developed a functional ex vivo RAD51-based test to identify HRD primary breast cancers. An integrated approach examining methylation, gene expression, and whole-exome sequencing was employed to ascertain the aetiology of HRD. Functional HRD breast cancers displayed genomic features of lack of competent HR, including large-scale state transitions and specific mutational signatures. Somatic and/or germline genetic alterations resulting in bi-allelic loss-of-function of HR genes underpinned functional HRD in 89% of cases, and were observed in only one of the 15 HR-proficient samples tested. These findings indicate the importance of a comprehensive genetic assessment of bi-allelic alterations in the HR pathway to deliver a precision medicine-based approach to select patients for therapies targeting tumour-specific DNA repair defects. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

The genetic landscape of breast carcinomas with neuroendocrine differentiation.

Neuroendocrine breast carcinomas (NBCs) account for 2-5% of all invasive breast cancers, and are histologically similar to neuroendocrine tumours from other sites. They typically express oestrogen receptor (ER), and are HER2-negative and of luminal 'intrinsic' subtype. Here, we sought to define the mutational profile of NBCs, and to investigate whether NBCs and common forms of luminal (ER+ /HER2- ) breast carcinoma show distinct repertoires of somatic mutations. Eighteen ER+ /HER2- NBCs, defined as harbouring >50% of tumour cells expressing chromogranin A and/or synaptophysin, and matched normal tissues were microdissected and subjected to massively parallel sequencing targeting all exons of 254 genes most frequently mutated in breast carcinomas and/or related to DNA repair. Their mutational repertoire was compared with that of ER+ /HER2- breast carcinomas (n = 240), PAM50-defined luminal breast carcinomas (luminal A, n = 209; luminal B, n = 111) and invasive lobular carcinomas (n = 127) from The Cancer Genome Atlas. NBCs were found to harbour a median of 4.5 (range 1-11) somatic mutations, similar to that of luminal B breast carcinomas (median = 3, range 0-17) but significantly higher than that of luminal A breast carcinomas (median = 3, range 0-18, p = 0.02). The most frequently mutated genes were GATA3, FOXA1, TBX3, and ARID1A (3/18, 17%), and PIK3CA, AKT1, and CDH1 (2/18, 11%). NBCs less frequently harboured PIK3CA mutations than common forms of ER+ /HER2- , luminal A and invasive lobular carcinomas (p < 0.05), and showed a significantly higher frequency of somatic mutations affecting ARID1A (17% versus 2%, p < 0.05) and the transcription factor-encoding genes FOXA1 (17% versus 2%, p = 0.01) and TBX3 (17% versus 3%, p < 0.05) than common-type ER+ /HER2- breast carcinomas. No TP53 somatic mutations were detected in NBCs. As compared with common forms of luminal breast carcinomas, NBCs show a distinctive repertoire of somatic mutations featuring lower frequencies of TP53 and PIK3CA mutations, enrichment for FOXA1 and TBX3 mutations, and, akin to neuroendocrine tumours from other sites, ARID1A mutations. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Genetic analysis of microglandular adenosis and acinic cell carcinomas of the breast provides evidence for the existence of a low-grade triple-negative breast neoplasia family.

Acinic cell carcinoma is an indolent form of invasive breast cancer, whereas microglandular adenosis has been shown to be a neoplastic proliferation. Both entities display a triple-negative phenotype, and may give rise to and display somatic genomic alterations typical of high-grade triple-negative breast cancers. Here we report on a comparison of previously published data on eight carcinoma-associated microglandular adenosis and eight acinic cell carcinomas subjected to targeted massively parallel sequencing targeting all exons of 236 genes recurrently mutated in breast cancer and/or DNA repair-related. Somatic mutations, insertions/ deletions, and copy number alterations were detected using state-of-the-art bioinformatic algorithms. All cases were of triple-negative phenotype. A median of 4.5 (1-13) and 4.0 (1-7) non-synonymous somatic mutations per carcinoma-associated microglandular adenosis and acinic cell carcinoma were identified, respectively. TP53 was the sole highly recurrently mutated gene (75% in microglandular adenosis versus 88% in acinic cell carcinomas), and TP53 mutations were consistently coupled with loss of heterozygosity of the wild-type allele. Additional somatic mutations shared by both groups included those in BRCA1, PIK3CA, and INPP4B. Recurrent (n=2) somatic mutations restricted to microglandular adenosis or acinic cell carcinomas included those affecting PTEN and MED12 or ERBB4, respectively. No significant differences in the repertoire of somatic mutations were detected between microglandular adenosis and acinic cell carcinomas, and between this group of lesions and 77 triple-negative carcinomas from The Cancer Genome Atlas. Microglandular adenosis and acinic cell carcinomas, however, were genetically distinct from estrogen receptor-positive and/or HER2-positive breast cancers from The Cancer Genome Atlas. Our findings support the contention that microglandular adenosis and acinic cell carcinoma are part of the same spectrum of lesions harboring frequent TP53 somatic mutations, and likely represent low-grade forms of triple-negative disease with no/minimal metastatic potential, of which a subset has the potential to progress to high-grade triple-negative breast cancer.

HMGA1 Expression in Human Hepatocellular Carcinoma Correlates with Poor Prognosis and Promotes Tumor Growth and Migration in in vitro Models.

BACKGROUND: HMGA1 is a non-histone nuclear protein that regulates cellular proliferation, invasion and apoptosis and is overexpressed in many carcinomas. In this study we sought to explore the expression of HMGA1 in HCCs and cirrhotic tissues, and its effect in in vitro models. METHODS: We evaluated HMGA1 expression using gene expression microarrays (59 HCCs, of which 37 were matched with their corresponding cirrhotic tissue and 5 normal liver donors) and tissue microarray (192 HCCs, 108 cirrhotic tissues and 79 normal liver samples). HMGA1 expression was correlated with clinicopathologic features and patient outcome. Four liver cancer cell lines with stable induced or knockdown expression of HMGA1 were characterized using in vitro assays, including proliferation, migration and anchorage-independent growth. RESULTS: HMGA1 expression increased monotonically from normal liver tissues to cirrhotic tissue to HCC (P<.01) and was associated with Edmondson grade (P<.01). Overall, 51% and 42% of HCCs and cirrhotic tissues expressed HMGA1, respectively. Patients with HMGA1-positive HCCs had earlier disease progression and worse overall survival. Forced expression of HMGA1 in liver cancer models resulted in increased cell growth and migration, and vice versa. Soft agar assay showed that forced expression of HMGA1 led to increased foci formation, suggesting an oncogenic role of HMGA1 in hepatocarcinogenesis. CONCLUSIONS: HMGA1 is frequently expressed in cirrhotic tissues and HCCs and its expression is associated with high Edmondson grade and worse prognosis in HCC. Our results suggest that HMGA1 may act as oncogenic driver of progression, implicating it in tumor growth and migration potential in liver carcinogenesis.

Massively parallel sequencing of phyllodes tumours of the breast reveals actionable mutations, and TERT promoter hotspot mutations and TERT gene amplification as likely drivers of progression.

Phyllodes tumours (PTs) are breast fibroepithelial lesions that are graded based on histological criteria as benign, borderline or malignant. PTs may recur locally. Borderline PTs and malignant PTs may metastasize to distant sites. Breast fibroepithelial lesions, including PTs and fibroadenomas, are characterized by recurrent MED12 exon 2 somatic mutations. We sought to define the repertoire of somatic genetic alterations in PTs and whether these may assist in the differential diagnosis of these lesions. We collected 100 fibroadenomas, 40 benign PTs, 14 borderline PTs and 22 malignant PTs; six, six and 13 benign, borderline and malignant PTs, respectively, and their matched normal tissue, were subjected to targeted massively parallel sequencing (MPS) using the MSK-IMPACT sequencing assay. Recurrent MED12 mutations were found in 56% of PTs; in addition, mutations affecting cancer genes (eg TP53, RB1, SETD2 and EGFR) were exclusively detected in borderline and malignant PTs. We found a novel recurrent clonal hotspot mutation in the TERT promoter (-124 C>T) in 52% and TERT gene amplification in 4% of PTs. Laser capture microdissection revealed that these mutations were restricted to the mesenchymal component of PTs. Sequencing analysis of the entire cohort revealed that the frequency of TERT alterations increased from benign (18%) to borderline (57%) and to malignant PTs (68%; p < 0.01), and TERT alterations were associated with increased levels of TERT mRNA (p < 0.001). No TERT alterations were observed in fibroadenomas. An analysis of TERT promoter sequencing and gene amplification distinguished PTs from fibroadenomas with a sensitivity and a positive predictive value of 100% (CI 95.38-100%) and 100% (CI 85.86-100%), respectively, and a sensitivity and a negative predictive value of 39% (CI 28.65-51.36%) and 68% (CI 60.21-75.78%), respectively. Our results suggest that TERT alterations may drive the progression of PTs, and may assist in the differential diagnosis between PTs and fibroadenomas. Copyright © 2015 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Massively parallel sequencing analysis of synchronous fibroepithelial lesions supports the concept of progression from fibroadenoma to phyllodes tumor.

Phyllodes tumors (PTs) and fibroadenomas (FAs) are fibroepithelial lesions (FELs) of the breast. Although mutations affecting exon 2 of MED12 are highly recurrent in FAs and PTs, TERT promoter hotspot mutations are frequently found in PTs but are vanishingly rare in FAs. Malignant transformation of benign PTs is well-documented, but the progression from FA to PT remains a matter of contention. Here we report on the somatic genetic alterations in multiple ipsilateral synchronous FELs (three FAs, one benign PT, and one malignant PT) occurring in the same patient. DNA samples extracted from each tumor and matched normal tissue were subjected to targeted massively parallel sequencing using the Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets (MSK-IMPACT) assay. This analysis revealed MED12 mutations in all lesions. One FA and the benign PT harbored a MED12Gly44Val mutation, whereas another FA and the malignant PT displayed a MED12Gly44Asp mutation. The remaining FA had an independent distinct MED12Gly44Cys mutation. A formal clonality analysis suggested a clonal relationship between the FELs with identical MED12 mutations (P<0.05). A clonal TERT promoter hotspot mutation was identified exclusively in the malignant PT. The identification of distinct MED12 mutations in multifocal ipsilateral and synchronous FELs supports the notion that co-existing mammary fibroepithelial tumors can arise independently. Conversely, the co-existence of identical MED12 mutations indicates clonal relatedness among FAs and PTs, corroborating the hypothesis that FAs may constitute the substrate from which PTs develop. Our findings also support the notion that acquisition of TERT promoter mutations may drive the progression of FELs.

High-resolution genomic profiling of thyroid lesions uncovers preferential copy number gains affecting mitochondrial biogenesis loci in the oncocytic variants.

Oncocytic change is the result of aberrant mitochondrial hyperplasia, which may occur in both neoplastic and non-neoplastic cells and is not infrequent in the thyroid. Despite being a well-characterized histologic phenotype, the molecular causes underlying such a distinctive cellular change are poorly understood. To identify potential genetic causes for the oncocytic phenotype in thyroid, we analyzed copy number alterations in a set of oncocytic (n=21) and non-oncocytic (n=20) thyroid lesions by high-resolution microarray-based comparative genomic hybridization (aCGH). Each group comprised lesions of diverse histologic types, including hyperplastic nodules, adenomas and carcinomas. Unsupervised hierarchical clustering of categorical aCGH data resulted in two distinct branches, one of which was significantly enriched for samples with the oncocytic phenotype, regardless of histologic type. Analysis of aCGH events showed that the oncocytic group harbored a significantly higher number of genes involved in copy number gains, when compared to that of conventional thyroid lesions. Functional annotation demonstrated an enrichment for copy number gains that affect genes encoding activators of mitochondrial biogenesis in oncocytic cases but not in their non-oncocytic counterparts. Taken together, our data suggest that genomic alterations may represent additional/alternative mechanisms underlying the development of the oncocytic phenotype in the thyroid.

SF3B1 mutations constitute a novel therapeutic target in breast cancer.

Mutations in genes encoding proteins involved in RNA splicing have been found to occur at relatively high frequencies in several tumour types including myelodysplastic syndromes, chronic lymphocytic leukaemia, uveal melanoma, and pancreatic cancer, and at lower frequencies in breast cancer. To investigate whether dysfunction in RNA splicing is implicated in the pathogenesis of breast cancer, we performed a re-analysis of published exome and whole genome sequencing data. This analysis revealed that mutations in spliceosomal component genes occurred in 5.6% of unselected breast cancers, including hotspot mutations in the SF3B1 gene, which were found in 1.8% of unselected breast cancers. SF3B1 mutations were significantly associated with ER-positive disease, AKT1 mutations, and distinct copy number alterations. Additional profiling of hotspot mutations in a panel of special histological subtypes of breast cancer showed that 16% and 6% of papillary and mucinous carcinomas of the breast harboured the SF3B1 K700E mutation. RNA sequencing identified differentially spliced events expressed in tumours with SF3B1 mutations including the protein coding genes TMEM14C, RPL31, DYNL11, UQCC, and ABCC5, and the long non-coding RNA CRNDE. Moreover, SF3B1 mutant cell lines were found to be sensitive to the SF3b complex inhibitor spliceostatin A and treatment resulted in perturbation of the splicing signature. Albeit rare, SF3B1 mutations result in alternative splicing events, and may constitute drivers and a novel therapeutic target in a subset of breast cancers.

Benchmarking mutation effect prediction algorithms using functionally validated cancer-related missense mutations.

BACKGROUND: Massively parallel sequencing studies have led to the identification of a large number of mutations present in a minority of cancers of a given site. Hence, methods to identify the likely pathogenic mutations that are worth exploring experimentally and clinically are required. We sought to compare the performance of 15 mutation effect prediction algorithms and their agreement. As a hypothesis-generating aim, we sought to define whether combinations of prediction algorithms would improve the functional effect predictions of specific mutations. RESULTS: Literature and database mining of single nucleotide variants (SNVs) affecting 15 cancer genes was performed to identify mutations supported by functional evidence or hereditary disease association to be classified either as non-neutral (n = 849) or neutral (n = 140) with respect to their impact on protein function. These SNVs were employed to test the performance of 15 mutation effect prediction algorithms. The accuracy of the prediction algorithms varies considerably. Although all algorithms perform consistently well in terms of positive predictive value, their negative predictive value varies substantially. Cancer-specific mutation effect predictors display no-to-almost perfect agreement in their predictions of these SNVs, whereas the non-cancer-specific predictors showed no-to-moderate agreement. Combinations of predictors modestly improve accuracy and significantly improve negative predictive values. CONCLUSIONS: The information provided by mutation effect predictors is not equivalent. No algorithm is able to predict sufficiently accurately SNVs that should be taken forward for experimental or clinical testing. Combining algorithms aggregates orthogonal information and may result in improvements in the negative predictive value of mutation effect predictions.